Bohmian mechanics to high-order harmonic generation 
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Abstract 

This paper introduces Bohmian mechanics (BM) into the intense laser-atom physics to study 
high-order harmonic generation. In BM, the trajectories of atomic electrons in intense laser field 
can be obtained with the Bohm-Newton equation. The power spectrum with the trajectory of an 
atomic electron is calculated, which is found to be irregular. Next, the power spectrum associated 
with an atom ensemble from BM is considered, where the power spectrum becomes regular and 
consistent with that from quantum mechanics. Finally, the reason of the generation of the irregular 
spectrum is discussed. 
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I. INTRODUCTION 



Atom interacting with the intense laser field (ILF) can absorb more than one photon 
and then emits photons at harmonics of the laser frequency, which is well known as high- 
order harmonic generation (HHG) The harmonic spectrum has the following 
characteristics as a function of increasing frequency: a rapid decrease in harmonic intensity 
at first, followed by a plateau region with the harmonics having similar intensities; and then 
a rapid drop in harmonic strength beyond the plateau region. Recently, lots of theoretical 
works have studied this multiphoton phenomenon, including solving the time-dependent 
Schrodinger equation 

mm 

, the semiclassical trajectory methods, such as the three- 
step model 9| and the Feynman's path-integral approach in the strong field approximation 



step m 

10, lu 



12|, and even the classical trajectory method 13). 



Bohmian mechanics (BM) mm llq . or called quantum trajectory method }17| . is an- 
other alternative formalism of quantum mechanics. In the early years, BM has been used to 



re scattering by a square potential 
Recently, it has been applied to 



study some fundamental quantum phenomena such as t 
barrier [18|, [19( and the diff 
study and analyze many different processes and phenomena in physics and chemistry, such 
as atom-surface scattering [21 



sociation of NOC1 and NO2 



221 ]. electron transport in mesoscqpic systems 



241 ]. and the chemical reactions 25 1. Also, BM in terms of 



quantum trajectory has been considered to study chaos 



26 



27 



23J, photodis- 



2a, |2jj. Recently, it has 



been successfully applied to the intense laser-atom physics to study the dynamics of the 
above-threshold ionization (ATI) photoelectron by authors j3o| . 

In this paper, BM is introduced to the intense laser-atom physics to study HHG, which 
allows us to follow the time evolution of individual electron trajectory in atomic system 
first. We then calculate the power spectrum with the electron trajectory and find the cor- 
responding spectrum is irregular. Next, we consider the power spectrum associated with an 
atom ensemble. In the atom ensemble, harmonic generation is determined by the Maxwell 
equation, where the electronic polarization of the atom ensemble plays a key role for the 
harmonic generation. We will show the electronic polarization of the atom ensemble gained 
from BM is equivalent to that obtained from quantum mechanics. On this condition, the 
power spectrum from BM is regular and consistent with that obtained from quantum me- 
chanics. Finally, we will briefly discuss the reason why an individual electron trajectory 
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generates the irregular spectrum. 

This paper is organized as follow: We will briefly introduce quantum trajectory method 
first. Then we show the Hamiltonian for the hydrogen atom in ILF. In Sec. IV, we calculate 
the harmonic spectrum from a single electron trajectory, following by harmonic generation 
associated with an ensemble of atoms interacting with the ILF in Sec. V. Finally, we briefly 
discuss the reason of the generation of the irregular spectrum with an individual electron 
trajectory and then conclude. 



II. QUANTUM TRAJECTORY FORMALISM OF BOHMIAN MECHANICS 



Quantum trajectory method comes from the following transformation of the time- 

n n 

dependent Schrodinger equation |14j . [15| . Firstly a wave function can be written in the 
polar form ip(r,t) = R(r,t) exp(iS(r,t)/h). Secondly ip(r,t) is substituted into the time- 
dependent Schrodinger equation. Then the real part of the resulting equation is 



and the imaginary part has the form 

% + V(pv) = 0, (2) 

where p(r, t) = R 2 (r,t), v = V5(r, t)/m, and Q(r,t) = —^^j^- These two equations 
look like the classical Hamilton- Jacobi equation and the equation of continuity, respectively, 
except Eq. ([1]) has an extra term Q which is usually called quantum potential in BM. Thus 
Eq. ([1]) is called the quantum Hamilton- Jacobi equation governed by an external potential 
V and the quantum potential Q. The motion of particle is guided by the Bohm- Newton 
equation of motion, md 2 r/dt 2 = — V(V + Q), or, equivalently, by 

dr/dt = VS(r,t)/m. (3) 

In practice, we solve the time-dependent Schrodinger equation to obtain i[>(r,t), and then 
S(r, t). In this way, the quantum trajectory of the particle can be gained by integrating Eq. 
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III. NUMERICAL SOLUTION OF THE TIME-DEPENDENT SCHRODINGER 
EQUATION 



The Schrodinger equation for Hydrogen atom in ILF can be written as (atomic units are 
used throughout) 



Here Ho(r) is the field-free Hydrogen atom Hamiltonian and V(r,t) is the intense laser- 
atom interaction: H (r) = —\-^ + jnz ~ ~, ^( r )^) = — r ■ E(t) = —zE(t), where the 
laser field is the linearly polarized field (E||z) and E(t) is the laser field profile. Due to 
the linearly polarized laser field, magnetic quantum number m of Hydrogen atom is a good 
quantum number, so that the problem of solving the time-dependent Schrodinger here can be 
simplified into a two-dimensional problem (see the inset in Fig. 1). In this work the solution 
of the time-dependent Schrodinger equation ip(r,t) is obtained by the grid method and the 
second-order split-operator method, which has been detailedly introduced by Tong et al 3l| . 

f E sin 2 (^) sm(u t), 0<t< 10T; 
The laser field profile is E{t) = { 20T ~ ~ where T = 2ir/u , 

I E sm(cu t), t>WT, 

and E and u are the electric field amplitude and angular frequency, respectively. We choose 
E = 0.0292 a.u. (i.e., J = 3.0 X 10 13 W/cm 2 ), and u = 0.0587 a.u. (i.e., A = 775nm). The 
initial wave function ip(r, 0) of the system is in the ground state of the field- free Hydrogen 
atom. 
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FIG. 1: Power spectrum of the electron with the initial position ro = (1.0 a.u., vr/4). The inset 
shows the direction of the laser beam and the electron coordinate of atomic system. 
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IV. HARMONIC SPECTRUM ASSOCIATED WITH A SINGLE ELECTRON 
TRAJECTORY 



After obtaining the wavefunction ip(r,t), we can follow the time evolution of electron 
trajector y r (t) by integrating Eq. ([3]) with the electron initial position r . According to 
BM J14I. I15I. 16], the initial electron distribution is \ip(r, 0)| 2 . In this work, ip(r, 0) is the 
ground state of the field-free Hydrogen atom. Thus we can obtain an ensemble of electron 
trajectories with the corresponding electron initial positions. For an individual electron 
trajectory, the atomic dipole moment [i z {t) along the laser polarization direction is z ■ r(i). 
The corresponding power spectrum is gained by the Fourier transformation of the dipole 



moment 



13, 



1 



t f - 1 



fi g (t)e 



—iu)t 



% Jt 



(4) 



where U — and the propagation time tf is 35T in this work. 

Explicitly, we take one atomic electron as an example with the initial position r = (1.0 
a.u., 7r/4) in polar coordinates. The electron trajectory r(i) is calculated by integrating Eq. 
([2]). The corresponding power spectrum is shown in Fig. 1 from Eq. (J3J). The spectrum 
has even harmonics, which is dominated by a Rayleigh scattering component at the laser 
frequency to = ujq. We have obtained lots of electron trajectories with different initial position 
ro according to the initial electron distribution \ip(r, 0)| 2 . The corresponding power spectra 
have the similar character described above. In the following, however, we will show the even 
harmonics of the power spectrum are coherently removed in an ensemble of atoms and the 
corresponding power spectrum is consistent with that obtained from quantum mechanics. 



V. HARMONIC GENERATION FROM AN ENSEMBLE OF ATOMS INTER- 
ACTING WITH THE INTENSE LASER FIELD 



The 



33. 



34j, 



rarmonic generation of an atom ensemble can be gained with the Maxwell equation 



35, 



V 2 S 



ld 2 £ And 2 ^ 



c 2 dt 2 c 2 dt 2 



where <^(r', t) is the electronic polarization of the atom ensemble at the point r' and t) 
describes the electromagnetic field (involving the fundamental and the harmonic fields). 
Obviously, the electronic polarization &(r',t) plays a key role for the harmonic generation 
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from the ensemble of atoms. We will show that the electronic polarization &(r',t) gained 
from BM is, approximately, equivalent to that obtained from quantum mechanics, i.e., both 
of their harmonic spectra from the atom ensemble are consistent. 

We first assume the number of atoms is N in a small volume AV near the point r' in the 
atom ensemble ( AV <C d 3 , where d is the focal spot size of laser). The electronic polarization 

N 

&(r') is defined as the atomic dipole moment per unit volume: &(r') = Y^Pi/AV, where 

i 

Pi is the atomic dipole moment of each atom in the small volume AV^. From the view 
point of quantum mechanics, the dipole moment of each atom is p(t) = (ip(r,t)\z\ip(r,t)) 
in the z direction, where r is the electron coordinate of atomic system and if)(r,t) is the 
electron wavefunction. Thus the electronic polarization of the atom ensemble obtained from 
quantum mechanics is 

&>(T',t) = Np(t)/AV (5) 

in the z direction. 




FIG. 2: A small volume AV near the point r' in the atom ensemble, which includes lots of atoms 
with the different dipole moments from the view of BM. The inset shows the electron coordinate 
of atomic system. 



culate the electronic polarization of the atom ensemble from 



On the other hand, we ca 
BM. According to BM Q,[l5, [l6|, the probability that an electron lies between r and r + dr 
in the atomic system at the time t is given by \ij)(r, t)| 2 dr, where the corresponding atomic 
dipole moment is r (in atomic units). Then in the small volume AV of the atom ensemble, 
the number of atoms with atomic dipole moment r is (N\ip(r, t)| 2 dr). The sum of all kinds 



of atomic dipole moments in the volume AV is ^'lV'( r i) ^)\ 2( ^ r j r j- Thus the electronic 



j 



polarization £^(r',t) at the point r' in the atom ensemble is ^2 N\ip(rj,t)\ drjTj/ AV (see 



j 



Fig. 2). For an atom ensemble, we can replace the sum by an integral, 

E AV ~* AV / ^ t){ Tdr - (6) 

Further, because the laser field is the linearly polarized field (E||z), the time-dependent 
wavefunction ip(r,t) is symmetrical along z axis. Thus in the z direction, 

^ ' 'r, t) \ 2 rdr = — — f \ip(r, t)\ 2 r cosOdr 



Al/./ AV 

N 



j |^(r,t)| 2 ^r, (7) 



AV 

where the last term equals Eq. (jSJ). In this way, the electronic polarization £P(r',t) gained 
from BM is, approximately, equivalent to that obtained from quantum mechanics, i.e., both 
of their harmonic spectra from the atom ensemble are the same. Figure 3 shows the har- 
monic spectrum from an ensemble of Hydrogen atoms interacting with ILF (laser confocal 
parameter b = 5.0 mm, gas density 10 17 atoms/cm 3 described by a truncated Lorentzian in 
the x direction with L = 0.5 mm [33j), which has a clear plateau region and cutoff at the 
11th order. 




FIG. 3: Harmonic spectrum associated with an ensemble of Hydrogen atoms at the laser field 
I = 3.0 x 10 13 W/cm 2 and A = 775nm. 
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VI. DISCUSSION 



In BM, harmonic spectrum from a single electron trajectory has even harmonics (see 
Fig. 1), but they are removed in the spectrum associated with an ensemble of atoms (see 
Fig. 3). The reason is that an individual electron trajectory does not possess an inversion 
center but the electron trajectories of the atom ensemble have such centre. Here we 
will numerically show that two electron trajectories are enough to get rid of the unphysical 
even harmonics if the initial positions ro of the two electrons are centrally symmetric. Let's 
take two electrons as an example with the initial positions Tq = (1.0 a.u., vr/4) and rjj = (1.0 
a.u., 57r/4), respectively. The corresponding atomic dipole moment /i z (t) along the laser 
polarization direction is z ■ [r 1 (t) + r 2 (t)] and the power spectrum gained from Eq. (4) is 
shown in Fig. 4 (solid curve), which has only odd-order harmonics. In addition, we calculate 
the HHG power spectrum in the length form from the time-dependent Schrodinger equation 
(dotted curve in Fig. 4) [3l|. Note that these two curves can basically overlap if the value of 
the power spectrum from BM is multiplied by a factor of 0.4, i.e., both of them have the same 
relative intensity. This is an interesting result, but the reason why HHG power spectra from 
BM and the time-dependent Schrodinger equation have the same relative intensity should 
be studied in the future. 




6 8 10 12 14 16 18 

Harmonic order (Wm ) 



FIG. 4: The solid curve corresponds to power spectrum of two Bohmian electrons with the centrally 
symmetric initial positions: = (1.0 a.u., vr/4) and rjj = (1.0 a.u., 5tt/4). The dot is the HHG 
power spectrum in the length form from the time-dependent Schrodinger equation (TDSE). 
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VII. SUMMARY 



In summary, we introduce BM into the intense laser-atom physics to discuss HHG. It 
allows us to follow the time evolution of each electron trajectory in an atomic system. We 
find that the power spectrum from an individual electron trajectory has the even unphysical 
harmonics. But the even harmonics are coherently removed in the ensemble of atoms and 
the power spectrum is consistent with that obtained from quantum mechanics. The reason 
of the appearance of the even harmonics is an individual electron trajectory does not possess 
an inversion center, but the electron trajectories of the atom ensemble have such centre. 
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